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CONFIDENTIAL. 


R.  C.  279. 


MINISTRY  OF  HCME  SECURITY. 

CIVIL  DEFENCE  RESEARCH  COMMITTEE. 

Notea  on  H.A.  Bethel  "Theory  of  armor  penetration1*. 

I.  Static  penetration. 

By  Professor  G.I.  Taylor,  F.R.S. 


The  first  part  of  this  paper  describes  the  static  stresses  in  a  long 
cylindrical  hollow  cylinder  and  in  a  flat  sheet  when  a  concentric  hole  is 
opened  out  by  radial  pressure  applied  over  its  surface.  Within  a  certain 
radius  the  material  is  assumed  to  be  overstrained  and  to  flow  radially. 

Outside  this  radius  the  conditions  are  elastic.  Far  the  thick  cylinder, 
where  it  is  assumed  that  there  is  no  extension  parallel  to  the  axis  of 
symmetry,  the  problem  and  its  solution  are  identical  with  those  given  in  text 
books  of  gunnery  in  connection  with  the  autofrettage  of  guns  and  with  those 
which  have  been  used  in  designing  cylinders  for  high  pressure  work.  In  this 
case  the  type  of  the  strain  can  be  related  immediately  to  a  single  variable, 
namely  the  radial  displacement  which  is  a  function  of  one  independent 
variable,  the  radius,  and  one  parameter,  the  radial  displacement  of  the  inner 
surface. 1  This  consideration  remains  true  when,  as  in  the  case  considered 
by  Dr.  Bethe,  the  strains  in  the  inner  plastic  region  are  not  small. 

The  hole  in  a  thin  plate  is  more  interesting  and  more  difficult  to 
analyse  because  it  is  no  longer  possible  to  treat  the  strain  as  two  dimens¬ 
ional,  so  that  the  relationship  between  plastic  strain  and  stress  must  be 
considered.  It  is  usually  assumed  that  hydrostatic  pressure  merely 
compresses  a  plastic  material  without  altering  its  strength  to  resist  shear 
stresses.  For  this  reason  it  is  sometimes  convenient  in  comparing  various 
theories  of  plasticity  to  use  reduced  principal  stresses  =  07  -  p, 

<£  =  <ra  -  p,  <r3'  =  <r3  -  p,  where  3p  =  cr,  +  +  cr  ,  so  that  07'  ♦  a'  +  <rJ  =  0. 

Similarly  reduced  principal  strains  e,'  =  e(  -  e ,  =  e2  -  e  ,  e^  =  e3  -  e 

where  e,  +  ez  +  e3  =  3e  and  e  represents  the  volumetric  strain.  The 
plasticity  relations  are  concerned  firstly  with  the  man-imnm  values  which  the 
stresses  07  ,  o^’  ,  <r,  can  attain  before  plastic  flow  occurs  and,  secondly, 
with  the  dependence  of  ej  ,  e^  ,  ej  on  07 1  ,  o ,  cr3  •  These  two  kinds  of 
plasticity  condition  are  quite  unrelated  to  one  another.  Of  the  first  type 
two  alternative  hypotheses  are  mentioned  by  Dr.  Bethe,  namely  those  of  Mohr 
and  v.  Mises,  and  he  points  out  that  there  is  but  little  difference  between 
them. 


For  two  dimensional  problems,  where  if  the  compressibility  be 
neglected  we  may  take  e3  =  0,  the  second  type  of  plasticity  condition  does 
not  affect  the  distribution  of  stress  in  the  plane  to  which  the  displacements 
are  confined.  This  is  because  when  e'3  =0,  e,'  =  -  ,  so  that  only  one 

kind  of  strain  is  possible  when  the  directions  of  the  principal  strains  are 
assumed  to  coincide  with  those  of  principal  stresses. 

The  case  is  very  different  when  the  strain  is  not  two  dimensional. 
Here  it  is  necessary  to  choose  seme  arbitrary  law  or  to  use  experimental  data. 
The  problem  can  be  visualised  by  thinking  of  the  relationship  between  the 
stress  ellipsoid  and  the  strain  ellipsoid. 

The  following  points  may  be  noticed 

(l)  The  absolute  magnitude  of  the  stress  ellipsoid  is  related  to  the 

strength  criterion,  e.g.  the  Mohr  or  v.  Mises  criteria. 


1 


not  counting  as  parameters  the  compressibility  of  the  material  or 
the  yield  strength. 


(2)  The  absolute  magnitude  of  the  strain  ellipsoid  bears  no  relation¬ 
ship  to  the  stresses  if  the  plastic  body  is  assumed  to  possess  the 
property  that  flow  will  occur  when  the  yield  stress  is  reached. 

(3)  it  is  a  necessaiy  condition  of  isotropy  of  the  plastic  material 
that  the  directions  of  the  principal  axes  of  the  stress  and  strain 
ellipsoids  shall  coincide. 

(4)  Owing  to  the  fact  that  ej  +  ea  +  ea  =  0  and  <r,  +  cr4'  +  o^(  =  0,  it  is 

necessary  to  know  only  one  relationship  in  order  to  determine  the  ratios 
ej  :  :  63  when  the  ratio  of  any  pair  of  ,  <ra ,  crj  is  known.  This 

relationship  can  conveniently  be  defined  in  terms  of  two  non— dimensional 
variables  p  and  V  (Lode's  variables) 
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<*■.  -  <r3 


1,  V  = 


2  e*.„T.  .g3 
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where  <T,  >  o^>  <5^.  These  variables  are  chosen  for  convenience  so  that 
p  lies  between  -1  and  +1 .  It  seems  that  all  plastic  materials  must 
satisfy  the  relationship  e,  >  ea>  e3  when  aj  >  o^>  <rt,  so  that  V  also 
lies  between  -1  and  +1.®  The  observed  relationship  between  p  and  V 
for  mild  steel,  soft  iron  and  copper  is  given  in  a  paper  by  Taylor  and 
Quinney^ ,  and  for  copper,  iron  and  nickel  by  Lode^.  For  all  these 
metals  the  relationship  is  substantially  that  shown  in  Fig.  1,  which  also 
contains  Taylor  and  Quinney's  experimental  results.  This  experimental 
relationship  may  be  compared  with  that  which  exists  in  all  Newtonian 
viscous  fluids,  namely  p  =  V  .  It  is  found  experimentally  that  p  = V 
for  lead  at  atmospheric  temperature  and  for  glass  heated  till  it  is  Just 
soft  enough  to  flow.  In  developing  theories  of  plasticity  seme  workers 
have  assumed  the  relationship  ft  *V  •  It  seems  unlikely  that  the 
divergence  between  the  observed  relationship  and  the  assumed  p  =  V 
will  give  rise  to  much  error  in  calculating  stress  and  strain  distri¬ 
butions.  It  will  be  noticed  that  the  assumption  p  =  V  ,  though  used 
by  v.  Mises,  is  quite  unrelated  to  v.  Mises'  criterion  of  strength.  The 
relationship  p  =  could  equally  well  be  used  with  Mohr's  strength 
relationship,  namely  that  flow  begins  when  cr,  -  0"3  =  constant  =  T. 


Bethe's  stress-strain  assumption. 

Bethe  considers  two  regions  of  plastic  flow,  the  outer  one 
extending  inwards  frem  the  outer  limit  of  plastic  flow  r  =  r,  to  the  radius 
r  =  at  which  the  tangential  stress  ceases  to  be  a  tension.  In  this  region 
the  radial  stress  must  be  taken  as  ,  the  tangent ial  tension  as  crj  ,  and  cTl  » 
the  intermediate  stress  normal  to  the  sheet,  is  zero.  Between  r  =  r,  and 
r  -  r,,  therefore.  Lode's  variable  p  is  positive  but  <1.  At  r  =  r^ ,  p  =  1 
since  at  that  point  (Ta  =  O',  =  0.  In  this  region  Bethe's  strain  assumption 
(which  he  attributes  to  Mohr)  is  that  the  plastic  flow  is  limited  to  the  plane 
of  the  sheet,  no  thickening  occurring  (see  p. 9  of  Bethe's  report).  If  the 
strain  is  limited  to  the  plane  of  the  sheet  ea  =  0  and  if  the  effect  of 
compressibility  is  neglected  e3  =  -e, .  Thus  in  the  region  r(>r  >ra,  V  =  0. 
This  is  shown  in  Fig.  1  by  means  of  the  line  AB. 

Though  Bethe's  strain  assumption  is  very  far  frcm  what  is  observed 
in  experiments  in  which  plastic  strains  are  measured,  yet  this  does  not 
necessarily  detract  frcm  the  value  of  his  calculation  of  stress  distribution 
in  the  region  r,  >r >ra,  because  with  the  "ideal"  plastic  body,  which  begins 
to  flow  as  soon  as  the  stress  reaches  a  given  value  and  continues  flowing 
until  the  stress  is  reduced,  an  infinitesimal  plastic  strain  may  enable  the 
equilibrium  stress  distribution  to  be  attained.  In  other  words,  if  only  a 
small  thickening  of  the  sheet  does  occur  it  will  produce  only  a  negligible 
effect  on  the  stress  distribution. 


1  Phil.  Trans,  Roy.  Soc. ,  230,  1 931 . 

2 

"Versuche  fiber  den  Einfluss  der  mittleren  Hauptspannung  auf  das 
Fliessen  der  Metalle,  Eisen,  Kupfer  und  Nickel",  Z.  Physik, 
vol.36  (1926). 
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In  the  range  r,  >r  >  r^,  when  the  maximum  stress  difference  is  that 
between  the  two  principal  stresses  in  the  plane  of  the  sheet,  the  equation  of 
equilibrium  is  sufficient,  with  Mohr's  strength  condition  prescribing  a 
constant  difference  between  them,  to  determine  the  stress.  Inside  the  radius 
rz,  i.  e.  when  rA>  r  >  b  where  b  is  the  radius  of  the  hole,  the  tangential 
stress  cannot  remain  positive  (tensile).  Two  alternatives  remain  - 


(a) 

00 


becomes  negative  (i. e.  there  can  be  a  compressive 

tangential  stress)  or 

=  0. 


a 

Bethe  rejects  alternative  (a)  because  in  that  case  <Tf  would  be  the 
intermediate  principal  stress  and  by  his  strain  assumption  it  would  be 
necessary  that  no  strain  could  take  place  in  the  tangential  direction.  This 
would  preclude  any  radial  displacement.  He  is  left  with  (b)  as  the  only 
possible  alternative  consistent  with  his  strain  assumption,  namely  (T9  =  0. 

This  alternative,  however,  suffers  from  very  severe  disadvantages.  The  stress 
at  every  point  is  one  which  is  symmetrical  about  the  radial  direction,  i. e. 
the  stress  ellipsoid  at  any  point  is  a  spheroid  whose  axis  of  symmetry  is 
along  a  radius.  On  the  other  hand  the  plastic  strain  which  according  to 
Bethe' s  calculation  results  from  this  symmetrical  or  uni -directional  stress  is 
very  far  from  symmetrical  and  is  variable  along  the  radius.  Expressed  in 
terms  of  Lode' s  variables  the  stress  in  the  range  r^>  r  >  b  is  represented  by 
f*-  -  1  while  the  strain  is  indeterminate  and  covers  a  range  of  the  line  U  *  1 
in  Pig.  1 . 


Since  the  alternative  (a)  that  0^  becomes  a  compressive  stress  when 
r  <  ra  is  perfectly  possible  if  other  stress-strain  assumptions  are  used,  it 
will  be  seen  that  the  sole  reason  for  Bethe 's  conclusion  that  =  0  is  that 
he  assumes  that  when  a  stress  is  applied  in  one  direction  (e. g.  a  pure 
pressure  or  tension  unaccompanied  by  transverse  stresses)  the  strain  is 
completely  indeterminate.  A  round  bar,  for  instance,  when  stretched  in  an 
ordinary  testing  machine,  would,  if  it  obeyed  Bethe 's  stress-strain  law,  in 
general  acquire  an  elliptical  section  and  it  is  this  assumed  asymmetrical 
property  of  plastic  material  which  alone  is  responsible  for  Bethe' s  conclusion 
that  0^=0. 


It  would  seem  better  to  abandon  the  attempt  to  give  a  reasoned 
justification  of  the  assumption  that  =  0  when  ra  >  r  >  b  and  to  fall  back  on 
the  fact  that  this  assumption  enables  a  stress  distribution  to  be  determined 
without  reference  to  the  strain.  The  equilibrium  equation  then  suffices  to 
determine  the  thickness  of  the  plate.  Comparison  between  the  results 
obtained  by  assuming  that  0^=0  and  those  observed  experimentally  might  then 
afford  a  Justification  for  this  assumption  as  being  adequate  for  demonstrating 
the  features  of  the  mechanics  of  the  problem  which  do  not  depend  on  the 
relationship  between  plastic  stress  and  strain. 

Though  Bethe  manages,  by  endowing  his  plastic  material  with  the 
ability  to  suffer  unsymraetric  strains  when  subjected  to  a  symmetrical  stress, 
to  avoid  all  consideration  of  successive  steps  by  which  any  given 
configuration  of  finite  strain  is  attained,  this  simplification  cannot  in 
general  be  made.  In  fact,  so  far  as  I  am  aware,  no  problem  of  plastic  flow 
which  involves  finite  displacements  has  ever  been  obtained  except  in  cases 
such  as  the  expansion  of  an  infinite  cylindrical  tube  by  internal  pressure, 
where  synmetry  alone  enables  the  strain  to  be  determined.  For  this  reason  it 
seems  desirable  to  formulate  the  equations  for  plastic  radial  flow  round  a 
hole  in  a  sheet  in  a  farm  which  can  be  applied  to  any  desired  law  of  strength 
such  as  Mohr's  or  v.  Mises*  or  any  desired  relationship  between  Lode's 
variables  f  and  V  . 

Analysis  of  strain  round  an  expanding  radial  hole  in  a  sheet. 

When  a  hole  is  enlarged  the  finite  strain  at  apy  stage  is  made  up 
of  infinitesimal  elements  of  strain  which  vary  as  the  enlargement  proceeds. 

Thus  when  a  small  pin  hole  in  a  plate  is  enlarged  we  must  study  the  small 
strain  produced  in  an  element  of  the  sheet  which  was  originally  at  radius  s 
fran  the  pinhole,  when  the  hole  enlarges  from  radius  b  to  radius  b  +  $b. 
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In  the  more  general  case  when  the  Initial  radius  of  the  hole  in  the 
unstretched  sheet  is  not  zero  this  is  very  difficult  to  analyse,  but  when  the 
expansion  starts  from  a  small  pinhole  it  may  be  expected  that  the  configur¬ 
ation  when  the  hole  has  radius  b^  will  be  similar  to  that  round  the  hole  when 
its  radius  is  b,  except  that  the  radii  where  any  given  thickness  occurs  will 
be  changed  in* the  ratio  ba/b,.  Thus  if  h  is  the  thickness  and  u  the  radial 
displacement,  it  may  be  assumed  that  h/ho  and  u/b  and  also  the  stresses  are 
functions  of  s/b  only:  where  hc  is  the  initial  thickness  of  the  sheet. 

To  simplify  matters  I  have  assumed  that  the  compressibility  is  so 
that  it  may  be  neglected  and  the  material  taken  as  incompressible. 

The  relationship  between  the  small  strain  which  occurs  at  any  radius  during 
the  expansion  of  the  hole  through  a  small  increase  in  radius  from  b  to 
fc  +  §  b  can  be  understood  by  referring  to  Fig. 2.  Here  the  ordinates 
represent  u  and  the  abscissae  r. 


The  initial  radial  distance  s  of  the  element  which  at  a  subsequent 
stage  in  the  opening  out  of  the  hole  is  at  radius  r  is  related  to  u  by  the 
equation 

r  =  s  ♦  u.  .  0) 


In  Fig. 2,  therefore,  the  displacement  of  a  particle  from  its  initial  radius  s 
i*  represented  by  a  line  drawn  at  45°  to  the  axes.  In  particular  the 
displacement  of  the  particles  which  were  initially  at  the  pinpoint  where  the 
hole  began  is  represented  by  the  45  line  OPoP,  .  The  curved  line  P0AQe 
represents  the  relationship  between  r  and  u  which  it  is  the  object  of  the 
analysis  to  calculate.  At  a  subsequent  stage  of  the  expansion,  when  the 
hole  has  expanded  from  radius  b  to  radius  b  +  $  b,  the  curve  P,  B  C  Q , 
representing  displacement  is  similar  to  P„  AQ0  tut  with  its  linear  dimensions 
increased  in  the  ratio  (b  +  $b)  :  b;  thus  in  Fig.  2  fi  Po  _  AC  _  AD  _  J_b 
so  that  0Po  AO  r  h 

AD  =  r§t/t°  .....  (2) 


If  Sr  is  the  change  in  r  for  a  given  particle  of  material  when  the  Qhole 
expands  from  b  to  b  +  §  b,  $  r  is  found  by  drawing  the  line  AB  at  45  to  the 
axes  to  meet  the  curve  P,  BCQj  in  B.  If  Sb/b  is  small  enough,  the  arc  CB 
may  be  taken  as  straight  so  that  if  TT  -  <£  is  the  slope  of  CB  to  the  axis 


yp  =  -  tan  cC  .  (3) 

If  /8  is  the  angle  A  0  Q0,  tan  /3  =  u/r.  Frcm  the  gecmetiy  of  the  figure 
ABCD  (Fig. 2) 

Sr  =  AF=HF  =  CE  tan  ct  +  DA  tan/S  =  (DA-Sr)tanot  +DA  tan/3 

.  (4) 

Hence  .  . 

Sr  =  ftan*  ♦■*S%1)da  .  (5) 

\  1  +  tan  d  / 

and  frcm  (2) 

(n  \ 

Sr  =  J  r  &  (6) 

The  radial  strain  component  during  the  expansion  of  the  hole  frcm  b 
to  b  +  Sb  is  \  and  differentiating  (6)  with  respect  to  r  keeping 

$b  constant,  drv*r/ 
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Since  the  strain  during  expansion  of  the  hole  fran  b  to  b  ♦  Sb  is 
proportional  to  §b/b,  it  is  convenient  to  define  strain  components  £r,  c# 
and  Ez.  so  that  strains  during  the  small  enlargement  Sb  acre  $b/b, 
£^Sb/b,  ^Sb/b.  With  this  definition 


The  tangential  strain  is  simply 


£ 


9 


b  $r 
TP  r 


u  au 
r  ~  *r 

1  -  i!£ 

Sr  _ 


(8) 


(9) 


and  the  strain  perpendicular  to  the  sheet  is 


(10) 


The  thickness  h  at  any  stage  can  be  found  simply  frcm  the  equation  of 
continuity:  it  is  given  by 


_h 

ho 


O-sK1-^ 


....  (ii) 


where  hQ  is  the  initial  thickness  of  the  sheet. 

It  is  a  simple  matter  to  verify  that  (10)  is  consistent  with  (11). 

These  expressions  for  strain  take  simple  forms  when  expressed  in 
terms  of  a  new  independent  variable  ^  =  rx  and  a  new  dependent  variable 
tj  s  bx  *  (r-u)1.  Making  these  transformations  and  writing 


P 


(8)  and  (9)  became 


♦  iia 
vx 


....  (12) 

....  (13) 


S‘  =  1  ‘Ip 


(14) 


while  ( 1 1 )  reduces  to  the  simple  form  iv/h0  =  p.  ....(15) 

It  is  a  simple  matter  to  deduce  (lif)  directly  from  (15). 

The  stress  equilibrium  equation  for  a  thin  sheet  is 


h (<Ty  ~  (Tj) 

r 


(16) 
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Two  possible  alternative  forms  for  the  strength  condition  might 
be  considered: - 


(a)  Mohr's  stress  criterion  which  may  be  written 

CT0  -  =  T  if  0^  is  positive,  i.e.  tensile' 

or  -  cr  *  T  if  <f9  is  negative,  i.e. 


(17) 


compressive 


(b)  Mises'  condition  which  may  be  written,  when  <TX  =  0, 

♦  0 (Ty.  =  constant  •••»  0®) 

This  reduces  to  -  OV  =  constant  if  CTfi  =0  and  so  is  identical  with 
Mohr's  in  that  case. 

If  Bathe 1 s  assumption  that  0^  =  0  combined  with  (Tr  *  constant  is 
used,  (16)  leads  to 


hr  »  constant  =  . ...  (19) 

where  rjt  is  the  outer  boundary  of  the  region  of  finite  plastic  strain. 
Substituting  in  ( 1 1 ) 

?  -0-s)0-fc) 

which  gives  on  integration 

£(r-u)X  »  rrx  +  constant 

Since  u  =  0  when  r  =  ra  the  constant  is  -^{ri)iand 

U  =  r  -  j(2r  -  ra)ra 

The  inner  boundary  is  where  b  =  r  =  u,  so  that  from  (22) 

b  =  ^ra  ....  (23) 

which  is  Bathe's  result  if  ra  is  identified  with  his  rA. 


....  (20) 

....  (21) 

....  (22) 


Plastic  strain  assumption  ^  =  V. 

The  only  simple  law  so  far  proposed  for  the  relationship  between 
plastic  stress  and  strain  which  is  consistent  with  isotropy  and  at  the  same 
time  resembles  what  is  observed  with  metals  is  that  represented  in  Lode's 
variables  by  \A-  a  V ,  i.  e.  small  strains  or  rates  of  extension  in  the 
principal  directions  are  proportional  to  OJ  -  p,  (Tz  -  p,  0^  -  p.  In  the 
present  case  where  <3^=0  this  is  represented  by  the  equation 


0> 

O'r 


£»-£■ 

TXT-. 


+  2£e 

2Er  *  t9 


(24) 


Substituting  (15)  and  (24)  in  (16)  the  equilibrium  condition  reduces  to 

2  Is  (pa.)  *  -  0  (25) 
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This  equation  must  be  used  in  conjunction  with  a  strength  criterion.  Mohr's 
criterion  (a)  will  be  used.  In  this  case  (25)  assumes  two  different  forms 
according  as  is  negative  (i. e.  compressive)  or  positive  (i. e.  tangential 
tension).  These  are 


<f9  negative,  OJ.  =  -Y  so  that  (25)  becomes 

211  *  f  (■&=&■) 

C "  positive,  CT*.  -  <S~9  a  -Y  so  that  from  (24)  <*>. 


(26) 


(-r) 


....  (27) 


Substituting  for  Sr  and  from  (13)  and  (14)  the  resulting 
equations  may  be  written: - 

C Tp  negative  (tangential  ccmpresaion) 


,x(p?)  *  q(-’  *%)  *§)  =  0 


In  this  case 


if . 

-2q  +  ^  * 


£2  _  -  .  la! 

<ry  "  P 


and  in  terms  of  Mohr' s  strength  criterion  the  stresses  are 

*  *  -r’  *  =  -t(£) 

CTfl  positive  (tangential  tension) 


...  (28) 


...  (29) 


....  (30) 


$P  J 


+ 


£-2a  +  -?i  =  0 

*  s1  *jp 


(3D 


where  w  is  written  for  dq 

dj 


The  expressions  for  £©/€».  and  0^/ (TV  cannot  be  simplified  by 
using  the  equation  of  equilibrium  and  the  full  expressions  derived  from. 
(13),  04)  and  (24)  must  be  used,  namely 


e  1  -2L 

L*  .  yp _ Sk 

t*  _4  7  233,  ~  JL  °r 

1  +  p*  +  5p 


- £22 -  -  1  _  (32) 

4^q  ♦  ^  -  P* 


and  in  terms  of  Mohr's  condition  -  0"3  a  Y 


the  stresses  are  now 
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<rr  =  -Y/(l  -  <r,/a;)  ) 

Oi  -  ^)/(i  -«/«;>  1  <33) 

It  will  be  seen  that  (31)  is  an  ordinary  differential  equation  of 
the  third  order  and  first  degree  while  (28)  is  of  the  second  order  and 
second  degree.  The  reason  for  this  difference  lies  in  the  form  of  Mohr's 
strength  condition.  When  a>  is  positive  three  boundary  conditions  can  be 
assigned  at  any  given  value  (i. e.  of  r).  These  might,  for  instance;,  be 
u/r,  h/he  and  which  can  be  transformed  directly  in  assigned  values  of  q, 

p  and  rt  .  'When  <Te  is  negative  (Tt  cannot  be  assigned  arbitrarily;  it  is 
in  fact  constant.  Thus  only  p  and  y  can  be  assigned  arbitrarily. 


Boundary  condition  at  the  elastic-plastic  boundary. 

The  elastic  stresses  due  to  radial  displacement  in  an  infinite 

sheet  are 

-<*  =  (Tg  =  t*/tx  ....  (34) 

where  r,  is  the  radius  at  which  <£  -  <5^  =  -Y.  The  corresponding  small 
radial  displacement  is 

....  (35) 


where  E  is  Young's  Modulus  and  m  is  Poisson's  ratio.  In  the  present 
investigation  compressibility  will  be  neglected  and  we  will  take  m  = 

In  the  elastic  region  therefore 

V  *  r‘(l  -f)‘  .  3  (l  -Ilf1)  -•  (36) 

At  the  inner  boundary  of  the  elastic  region  therefore 

p  =  dy/d'S  si,  q  =  o,  ....  (37) 

At  the  outer  boundary  of  the  plastic  region  since  Oe  is  positive  Mohr's 
criterion  ensures  that  0~r  -  CT0  =  -Y.  Since  OV  is  necessarily  continuous 
through  r  =  r,  ,  and  it  is  assumed  that  o;  -  *  -Y  at  the  elastic  limit 

in  the  elastic  region,  <Te  must  be  continuous  through  r  »  r,  and  equal  to 
^Y.  It  is  important  to  notice  the  reason  why  <Te  is  continuous  at  the 
plastic  boundary  in  this  case,  because  it  is  not  necessary  in  general  that 
a},  shall  be  continuous  when  Mohr's  criterion  is  used.  It  will  be  shown 
in  fact  that  (To  is  discontinuous  on  the  circle  r  =  rA  within  the  plastic 
region  where  (Te  *  0. 

Strains  and  displacements  when  r,  >  r  >  rx. 

In  the  region  within  the  circle  r  a  r(  where  is  positive,  it 
will  be  found  that  the  strains  are  small,  being  of  order  Y/E.  Assuming 


that  7]  =  ^  (1 

-a^()  and  p  =  1  +«p,  , 

where  «■  =  3Y/2E 

P  = 

fi-i  -  ff' 

so  that  p(  =  -  f),  j-2' 

....  (38) 

and  q  = 

....  (39) 

_ 


u 


m 
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When  o~e  is  positive 


♦  V3T  ♦  4^7q/p  ) 

+  V/^  +  Tjq/V  / 


....  (40) 


Substituting  from  (38)  and  (39)  in  (40)  and  neglecting  terms  in  otx  compared 
with  those  containing  <X 


g-.  .  -ft  (<>♦», 4V)  ....  (M) 

where  y  =  d^,/d>f  and  if'**  d.YVdT* 

Substituting  this  in  (27) 


0 


...*  (42) 


Neglecting  terms  which  contain  <X  as  a  factor  compared  with  those  that  do 
not,  p  may  be  taken  as  1  and  (42)  may  then  be  integrated  giving 


♦  KnS 


constant 


(43) 


The  boundary  conditions  at  r  =  r(  are  Pj  =  0,  Vt  s  0,  q  =  0,  so  that 
y'X  _  '$x(y'=  -2W'5=  2.  The  constant  in  (43)  is  therefore 

-1  +  i6i(r,x).  Writing  3  for  -£n(y/r,A)»  (43)  therefore  becomes 

*  0  -•  (w0 

The  integral  of  (44)  is 

■tn  (f»  +  $  +  3#i(3  +5)  =  constant 

or  £  1  +  r,a)j  =  -1  ••••  (45) 

the  constant  being  chosen  so  that  the  correct  boundary  conditions  are 
satisfied  at  r  =  r, . 

Since  if'  =  dty/d'S  ,  (45)  may  130  integrated.  Hence 


V,  =  - - 


2(1 


....  (46) 


The  equation  for  the  thickness  is 

k  "  3s"  1  -a(vff)  =  1 

The  displacement  is 

|  y 

u  =  ~  rf*  =  ^17  r  a  i«r 


+  <x 


1  1+5 

7  ‘  2(1  +  $*■)* 


|(i  ♦  w’S)  -  i 


....  (47) 


....  (48) 
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The  stress  can  be  found  by  substituting  frcm  (44) 


6  +S  1  di|/  .  . 

-  3  ,  i.  e.  for -“^7“  in  (41 ) 

It  is  found  that 


<Tr  =  4Y(1  -S)  =  -2€n(r/r,)} 

....  (49) 

and  hence 

O'*  =  r  -  <r„  =  *Y  { 1  +  2Mr/rf )} 

....  (50) 

This  is  the  well  known  result  which  can  be  obtained  without  considering  the 
strains  and  displacements,  if  it  is  assumed  that  the  thickness  of  the  plate 
does  not  vary* 

Values  when  crs  =  0. 

The  radius  r^  at  which  g^=  0  is  frcm  (50)  ra  =  r,/  -/e  =  0.6o6r( 
and  corresponds  with  5  =  -1 . 

Though  the  stress  distribution  in  the  range  r,  >  r  >  ra  is 
identical  with  that  found  in  Bethe's  investigation  the  displacements  and 
strains  are  not  the  same.  For  the  case  when  Poisson's  ratio  is  £  Bethe 
finds  the  displacement  when  O~0  =  0  is 


uz  s  ?<xr(a/ri  =  iara(0)  =  £<*^(2.718)  ....  (51) 

Putting  5  =  -1  in  (4B)  the  displacement  according  to  the  present  strain 
hypothesis  is 

uz  ~  e4?a  ~  ~  =  **  (52) 


The  displacement  is  in  fact  about  6  per  cent,  greater  than  that  calculated 
on  Bethe's  strain  hypothesis. 

Putting  ^5  ■  -1  in  (47)  the  value  of  h/h,,  at  r  =  r,  is  1  +  joC 
and  frcm  (39)  and  (44) 

q .  -*(2^  g  'j .  -**(2  -  ft*)  = 

hence  frcm  (45)  * 

q=  ITTYTH-TSSF 


so  that  when  ^  =  -1 , 


q  = 


27 

T5 


o L  • 


....  (55) 
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Boundary  values  at  Og  =  0. 

At  the  circle  r  =  rx  where  cr^  =  0,  p  and  nq  are  continuous.  Just 
inside  the  circle  therefore,  where  crp  is  negative 


and 


+  50C 


....  (54) 


ffhen  c r9  is  negative  q  is  determined  "by  (28)  when  7/'5  and  P  &ro 
given,  Substituting  in  (28)  from  (54)  the  values  of  q^  found  by  solving 
the  resulting  quadratic  equation  are  (neglecting  terms  in  oc*) 


q^  =  -  1  -  202  and  q  ^  •••• 

Neither  of  these  values  is  the  same  as  q  s  -27<*/l6,  the  value  Just  outside 
the  boundary,  so  that  q  is  not  continuous  at  r  *  r^ . 

Trying  first  the  value  q'f  -  +  31<*/8>  ^  and  ^  a re  calculated 
fraa  (13)  and  (14).  Their  values  are 

£„  -  ♦  and  £,  =  -  27^ 

so  that  from  (24) 

£  -  -3 

and  since  <rr  is  continuous  and  equal  to  -Y,  CTg,  would  be  positive.  This  is 
inconsistent  with  the  condition  that  <Ta  is  negative  or  zero  inside  r  =  rA, 
thus  the  solution  q*5  =  +  31 V8  must  be  rejected.  The  only  alternative, 

=  _  1  202^  must  therefore  be  correct.  It  will  be  noticed  that  this 

4  64 

involves  a  discontinuity,  not  only  in  q  but  in  £*,  and  consequently  in  • 

This  discontinuity  arises  from  the  form  of  Mohr's  criterion.  It  would  not 

occur  if  von  Mises '  criterion  had  been  used. 

Discontinuity  in  £►■  and  gg  . 

Substituting  q  =  -A,  -  t  p  =  1  ♦  »  ^?/!f  =  1 

in  (13)  and  (14),  it  is  found  that  8®  =  "  27<<  »  -  ex  and 

substituting  these  in  (24)  8 

iS?  s  A  yi°<- 

<T».  2  16 

ffhen  <X  is  small,  i.  e.  when  E/Y  is  small,  we  may  neglect  oi  and  take  as  the 
boundary  condition  at  r  =  r^  for  calculating  the  stresses  and  displacements 
when  r  <  ra  the  values 

p  =  1 ,  “Vj  -  1»  -  (56) 

and  the  stresses  are  <T>.  =  -Y,  0~9  =  ^Y. 

Thus  the  stress  cTj,  suddenly  changes  frcm  0  to  a  compressive  stress 
of  ^Y  at  the  radius  r  =  ra. 
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Calculation,  of  stress  and  strain  when  r<  rA. 


To  calculate  the  distribution  of  stress  and  plastic  strain  inside 
the  radius  r  =  rx,  (28)  must  be  solved  step  by  step.  This  could  be  done 
using  the  boundary  values  ^/r*  =  1  =  23*/8>  p  =  1  +  3£*/2j  q'J  =  “£  - 8la/ 1 6 
for  any  given  value  of  ot  ,  but  it  will  suffice  for  the  present  work  to  neglect 
the  small  terns  containing  ot  and  take  as  the  boundary  values  ^  /r*  a  1 , 
pal,  qr^  a  -  If  s  J  is  the  magnitude  of  a  small  step,  the  correspond¬ 

ing  changes  p  and  ?]  may  be  taken  as 


$  ^  «  p§  )a 

§p  =  qS'S 


(57) 


After  calculating  the  values  of  p  and  7j  at  the  end  of  each  step  these 
values  are  inserted  in  (28)  and  the  resulting  quadratic  for  q  is  solved, 
the  root  which  derives  by  continuous  variation  of  7}  and  p  from  qr^  =  -  i 
being  chosen  in  each  case. 


The  results  of  application  of  this  process  are  given  in  Table  1 
and  are  shown  graphically  in  Pig.  3°  Values  of  the  principal  variables 
"5  /t*  and  V/rx  are  £P-V'8n  *-n  cols.  1  and  2,  Table  1.  Values  of  p  and  ~qr^ 
are  given  in  cols.  3  and  4»  Using  these  values  of  p,  q  and  ^  values  of 

<Te/<rr  calculated  from  (29)  are  given  in  col.  8,  and  the  corresponding 
values  of  C7^/T  and  (r8/T  in  cols.  9  and  10.  It  will  be  seen  that  <T0  which 
begins  as  a  compressive  stress  equal  to  half  the  radial  stress  at  the  outer 
limit  of  the  region  of  finite  plastic  flow  rapidly  decreases  till  when 

0.35  it  becomes  zero, and  if  the  process  is  carried  further,  using  (29), 
cre  becomes  a  tension.  When  J/r *  =  0.30,  for  instance,  the  calculated  value 

of  op/ 6"r  is  -0,1 24-  For  values  of  'j/r*  less  than  0.  35,  therefore,  the 
alternative  form  (31 )  of  the  equilibrium  equation  must  be  used. 

Since  is  continuous  and  equal  to  -Y  at  'jT/rf  a  0.35  and  »  © 
when  is  Just  greater  than  0.35,  while  Mohres  criterion  ensures  that 

a;  -  cr8  =  -Y  when  cB  is  positive,  it  seems  that  CT£  =  0  when  /v£  is  Just 
less  than  0.35°  Since  both  CTr  and  cre  are  therefore  in  this  case  continuous 
through  the  radius  where  <re  changes  sign,  £>•  and  £9  are  also  continuous. 

Hence  from  (13)  q  is  continuous.  The  values  of  rj ,  p  and  q  at  ^/r*  =  0.35 
can  therefore  be  inserted  in  (31 )  and  the  value  of  w  =  d’^/d'J3  at 
J  /r^  =0.35  determined. 

The  changes  in  ,  p  and  q  during  the  first  step  1x1  the  new 
region  are  calculated  using  the  formulae 

S’?  =  pS?  *  £q(SX)2  *  iw(St)3 

Sp  «  qST  ♦  MSf f’  •••*  (58) 

§  q  =  w  §  J 

Values  of  tj/t*,  p,  -qr*  and  wra4  found  in  this  way  are  given  in  the  lower 
part  of  Table  1,  corresponding  with  0.35  >!?>  >  0.205°  Values  of  , 

(Ty/ Y  and  o^/Y  from  (32)  and  (33)  are  given  fn  cols.  8,  9,  10  of  Table  1. 

Conditions  at  edge  of  hole. 

It  will  be  seen  in  Table  1  that  as  /raa  decreases  to  0. 21 ,  “-qrj 
and  wr2v  are  rising  very  rapidly.  A  study  of  the  values  of  the  terms  in 
(31)  reveals  that  by  the  time  =  0.21  is  reached,  one  term  on  the  R.H.S. 

of  the  equation  and  one  on  the  L.H.  S.  are  larger  than  any  other  terms.  The 
limiting  form  of  the  equation  when  7j  is  small  is  in  fact  < 


3wt? 

P 


-  2q 


....  (59) 
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» 

This  equation  can  be  integrated  twice,  thus 


and 


-  q  =  A  7  ^ 
p2-  =  B  -  6A7'J3 


} 


(60) 


A  sind  B  being  the  two  constants  of  integration.  To  determine  the  constants 
the  values  of  *} /r£  =0.012,  p  =  2.326,  qr*  =  -23.22  can  be  used  at  =  0.21. 
The  resulting  values  of  A  and  B  are 

A  =  1.2171*^,  B  =  6.78. 

The  limiting  value  of  p  when  T]  =  0  is  therefore 


p  =  \/ 6. 78  =2.61  ....  (6l) 

This  is  the  limiting  value  of  h/hc  at  the  edge  of  the  hole  and  may  Vie 
compared  with  Bethe's  value  2.0.  It  is  not  very  different  from  the  value 
at  /r£  =  0. 21 ;  to  find  the  limiting  value  of  >  therefore  it  is 
sufficient  to  take  p  as  constant  and  equal  to  2.  61  in  the  interval  during 
which  7  decreases  from  0.012  to  0.  Thus  the  limiting  value  of  /r* 
corresponding  with  the  edge  of  the  hole  is 


^  ^  r“Z  =  0.21  -  =  0.205  ....  (62) 

The  ratio  radius  of  finite  plastic  deformation  is 
radius  of  hole 

b^  Tea®?  =  2,21  ••••  (63^ 

This  may  be  compared  with  Bethe’s  value  2.0. 

Substituting  the  approximate  limiting  forms  of  p  and  q  from  (40) 
in  (32)  the  limiting  form  for  0~e/crY  is 

LVo  a?  =  4-  86  97 J+3?.  75  “  1 

This  tends  to  the  value  -1  as  indicated  in  the  last  figure  of  col.  8  and  the 
corresponding  values  of  <TY  and  are  therefore  -0.  5Y  and  +0.51* 

It  will  be  noticed  that  the  stress  at  the  internal  boundary  could 
have  been  predicted  a  priori  if  it  had  been  possible  to  assume  that  h/h0 
is  finite  at  r  =  b,  because  clearly  the  total  amounts  of  strain  in  the 
tangential  and  radial  directions  are  both  infinite  at  a  hole  which  has  been 
enlarged  from  a  pinhole.  Thus  the  state  of  strain  at  the  hole  is  such 
that  symmetry  alone  must  ensure  that  OJ  is  exactly  half  way  between  ov  and 
CT9  .  Since  05.  =  0,  0^.=  -  CTe  .  Similar  considerations  can  be  used  to 
understand  why  the  stress  at  points  Just  inside  the  boundary  r  =  rx 
corresponds  with  (o^/cr^)  =  +0.  5»  f or  at  the  edge  of  the  region  of  finite 
plastic  displacement,  where  the  radial  displacement  is  zero,  =  0.  Thus 
and  cre  must  therefore  be  exactly  half  way  between  a\.  and  cr^  . 
Hence,  since  =  0,  =  £<5;. 

Expressions  in  teims  of  radius  of  hole. 

The  radial  variable  is  expressed  in  terms  of  the  radius  of  the 
plastic  region.  To  express  the  results  in  terms  of  b,  it  is  necessary  to 
tabulate  r/b  =  2.2l/Y/ra.  These  values  are  given  in  col.  6,  Table  1. 

The  displacements  u/b  =  2.21  (J^  -  Jr)  )r^'  are  tabulated  in  col.  7. 
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The  radial  displacement  is  shown  graphically  in  Pig.  3  which  may  be 
compared  with  the  diagrammatic  sketch,  Pig. 2. 

Campari a on  with  Bethe's  results. 

The  thickness  ratio  h/h0  s  p  is  shown  in  Pig. 4  and  Bethe's  values, 
namely  fy/ho  =  2b/r,  are  also  shown.  It  will  be  seen  that  the  main  differ¬ 
ences  are  that  the  present  calculation  shows  the  "crater"  extending  further 
radially  than  Bethe's  and  at  the  same  time  the  "crater"  is  much  steeper  close 
to  the  hole.  At  first  sight  it  might  be  thought  that  the  extra  thickness  at 
r  =  b  above  Bethe's  2h0  means  that  the  work  done  in  expanding  the  hole  is 
greater  according  to  the  present  calculations  than  in  Bethe's  calculations, 
but  this  is  the  reverse  of  the  truth  for  the  radial  stress  at  the  hole  is  only 
-£Y  instead  of  Bethe's  -Y.  In  fact  the  work  done  in  expanding  to  a  given 

radius  is  only  2.61/4  =  0.65  of  the  work  done  if  Bethe's  strain  assumption  is 
used. 

Pig.  5  shows  the  distribution  of  stress.  This  is  of  course  very 
different  fran  Bethe's,  the  most  striking  difference  being  that  the  present 
calculations  predict  a  state  of  tangential  tension  in  a  ring  which  extends  to 
3C0S  of  the  radius  of  the  hole  from  its  edge  and  a  tangential  compression  from 
that  point  to  the  edge  of  the  region  of  large  plastic  distortion.  In  the 
plastic  region  r(>r  >rA,  where  smell  strains  comparable  with  the  elastic  strains 
occur,  the  stress  is  as  calculated  by  Bethe,  i„e.  there  is  a  tangential 
tension.  In  this  connection  it  may  be  noticed  that  in  comparing  calculations 
of  this  kind  with  the  behaviour  of  real  materials  a  metal  which  experiences 
considerable  hardening  with  cold  work  might  give  results  differing  widely  from 
the  above  theory.  The  extra  hardness  of  the  material  near  the  hole  might  be 
expected  to  prevent  the  formation  of  the  thin  lip  shown  in  Pig. 4  which  the 
analysis  predicts  for  an  "ideal"  plastic  solid. 

TABLE  1. 


1 

2 

5 

4 

5 

— S - 

7 

8 

9 

10 

SA* 

V«f 

P 

~^rx 

w*a 

r 

b 

u 

b 

<r> 

ov 

cr* 

Y 

Y 

1.0 

1.0 

1.0 

0.25 

2.21 

0 

+0.50 

-1.0 

-0.50 

\ 

0.90 

0.899 

1.025 

0.305 

2.096 

.001 

-1.0 

-0.47 

) 

0.80 

0.795 

1.055 

0.381 

1.978 

.008 

+0.440 

-1.0 

—0. 44 

)  ^ 

0.75 

0. 741 

1.075 

0.431 

1.915 

.012 

+0.397 

-1.0 
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